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ABSTRACT 

It is widely believed that the maximum energy of synchrotron photons when electrons are 
accelerated in shocks via the Fermi process is about 50 MeV (in plasma comoving frame). 
We show that under certain conditions, which are expected to be realized in relativistic shocks 
of gamma-ray bursts, synchrotron photons of energy much larger than 50 MeV (comoving 
frame) can be produced. The requirement is that magnetic field should decay downstream of 
the shock front on a length scale that is small compared with the distance traveled by the 
highest energy electrons before they lose half their energy; photons of energy much larger 
than 50 MeV are produced close to the shock front whereas the highest Lorentz factor that 
electrons can attain is controlled by the much weaker field that occupies most of the volume 
of the shocked plasma. 
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1 INTRODUCTION 

High energy radiation from many astrophysical objects is believed 
to be generated in shock heated plasma where charged particles are 
accelerated by the Fermi mechanism, eg. Fermi (1949), Axford et 
al. (1977), Bell (1978), Blandford & Ostriker (1978), Achterberg 
et al. (2001), Sironi & Spitkovsky (2011). The maximum particle 
energy is set by the condition that radiative losses between accel- 
eration episodes (which is the time it takes to travel from one side 
of shock front to the other) should be smaller than the energy gain. 
This upper bound to electron energy means that the emergent syn- 
chrotron radiation falls off to zero above a certain cut-off frequency 
which is about 50 MeV in the plasma rest frame independent of the 
magnetic field strength in the shock heated fluid (de Jager & Hard- 
ing, 1992; Lyutikov, 2009; Piran & Nakar, 2010; the derivation is 
provided in the next section). 

Many astrophysical objects, such as supernova (SN) remnants, 
pulsars,AGNs and gamma-ray bursts (GRBs), emit photons of en- 
ergy larger than ~ lO^MeV, and this upper limit has been used to 
eliminate the electron synchrotron process in shock-heated plasma 
as the radiation mechanism for these objects. However, we show in 
this Letter that the ~ 50 MeV upper limit can be violated under 
certain astrophysically realistic conditions. 

The upper limit of 50 MeV for synchrotron radiation is arrived 
at by assuming a uniform magnetic field throughout shock heated 
plasma. We show that this upper limit can be considerably raised 
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when the magnetic field is not uniform. A particular case we ana- 
lyze in some detail is where the magnetic field decays with distance 
down-stream of the shock front, and find that for a range of param- 
eters there can be a significant synchrotron flux even at several GeV 
(in fluid comoving frame). This is described in section 2. The calcu- 
lations presented in section 2 can be easily generalized to consider 
a situation where magnetic field fluctuates randomly instead of de- 
caying with distance from the shock front. The application to a few 
astrophysical systems are presented in §3. 



2 MAXIMUM SYNCHROTRON FREQUENCY FOR 
SHOCK HEATED PLASMA 

The argument that the maximum photon energy for synchrotron 
emission for shock accelerated electrons is ~ 50 MeV is straight- 
forward and goes as follows. 

Electrons gain energy by a factor ~ 2 whenever they cross 
a relativistic shock front and are scattered back to the other side. 
The time it takes for a charge particle of Lorentz factor 7^ to make 
this trip from one side of the shock front to the other is of order 
of the Larmor time, II = mcye/ (qB); where B is the comoving 
frame magnetic field strength, m is the particle mass, and q is its 
charge. The energy lost to synchrotron radiation during this time 
is SE ~ IlcttB'^^Ic/Q'k ~ arB'y^mc^ /{G-rrq). Particle accel- 
eration ceases when SE ~ mc^7e/2, and therefore the maximum 
Lorentz factor a particle can attain is given by 
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where ctt ~ 87rg^/(3m^c*) is the Thomson cross-section. The 
synchrotron photon energy corresponding to 'ymax is 
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Thus, the maximum photon energy for electrons is ~ 50 MeV and 
for protons ~ 10^ Ge\{3. These numbers might be overestimated 
by a factor 5-10 due to uncertainty regarding the time it takes par- 
ticles to travel from one side of the shock front to the other (which 
could be an order of magnitude larger than the Larmor time as- 
sumed in these calculations, e.g. Achterberg et al. 2001, Lemoine 
& Revenu 2006, Sagi & Nakar 2012, Uchiyama et al. 2007). All 
calculations in this paper are affected by this uncertainty. However, 
the ratio of Umax for the case of a uniform magnetic field to that 
for an inhomogeneous field geometry, discussed below, should be 
fairly secure. 

2.1 Inhomogeneous magnetic field and the maximum 
synchrotron frequency 

Let us consider that the magnetic field decays with distance down- 
stream of the shock front as 



B{x)^B^ix/Lp)-'' + Bo 



(3) 



where Lp is the length scale over which the field decays, Bu, & 
Bo are the strongest and weakest magnetic field strengths in the 
shocked fluid, and x ^ Lp. If magnetic field is generated by the 
Weibel mechanism then Lp is of order the plasma length scale 
(Medvedev & Loeb, 1999), i.e. 

Lp = [mpFsC'/i'iTvn.q^)]^^^ = 2.2xl0^cm(r,/rie)'/'' (4) 

where rrip is proton mass, Ts is the Lorentz factor of the shock front 
wrt the unshocked fluid, and is the number density of electrons 
in the shocked fluid comoving fram^. 

Particles accelerated in a shock influence the generation of 
magnetic fields and can substantially increase the length scale for 
field decay (Medvedev et al. 2005; Medvedev & Zakutnyaya, 2009; 
Keshet et al. 2009). A larger coherence length magnetic field pro- 
duced by high energy particles - which have larger plasma scale ~ 
decays on a longer time scale. However, these fields are also likely 
to be weaker (Medvedev & Zakutnyaya, 2009), although their true 
strength is quite uncertain. We, therefore, use measurements of 
magnetic fields in GRB afterglows for guidance. The volume aver- 
aged magnetic field energy density in GRB external shocks, while 
they are relativistic, is found to be a factor ~ lO'* smaller than the 
equipartition value for a good fraction of bursts (Panaitescu & Ku- 
mar, 2001, 2002; Bjomsson et al. 2004; Wei et al. 2006; Rykoff 
et al. 2006; Rol et al. 2007; Chandra et al. 2008; Liang et al. 2008; 
Gao et al. 2009; Xu et al. 2009; Cenko et al. 2010; Rossi et al. 201 1 ; 
Kumar & Barniol Duran, 2010). This suggests that the field decays 
with distance downstream by a factor ^ 10^ since the strength 
close to the shock front, due to the Weibel instability, is of order 
of the equipartition value. The length scale for this decay, however, 
is uncertain. For the sake of concreteness we take Lp to be plasma 

^ Although proton synchrotron process can produce 7-rays with energy of 
~ 10^ GeV in plasma comoving frame, this is a very inefficient mechanism 
and not likely to play a significant role in GRBs. 

^ We note that for relativistic shocks there is little difference between pro- 
ton and electron plasma length scales as long as electrons and protons ai'e 
in rough equipartition which is found to be the case from the study of nu- 
merous gamma-ray burst afterglows, e.g. Panaitescu & Kumar (2002). 



scale for our calculations. Fortunately, the results presented in this 
paper don't depend on the precise length-scale for field decay as 
long as it is much smaller than the distance traveled by the highest 
energy electrons before they lose their energy to radiation. 

The Larmor radius of an electron, RLije) ~ meC^7e/ (qB), 
increases with distance from the shock front as the magnetic field 
gets weaker (eq.[3}, and an electron traveling down-stream is likely 
to be sent back upstream when Rl f~x (x is the distance from the 
shock-front). Therefore, an absolute upper limit on 7e is set by the 
requirement that ii_L is smaller than the width of the shock heated 
plasma. 

A more stringent upper limit on 7e is obtained by the re- 
quirement that energy lost by an electron to radiation while trav- 
eling from up-stream to down-stream (in between two consecutive 
episodes of energy gain) should be smaller than me7eC^/2. 

The energy loss rate for an electron of LF 7e due to syn- 
chrotron radiation while traveling down-stream of the shock front 
is 
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We are interested in the case where Bw 3> Bo . However, since Lp 
- the length scale over which the magnetic field decays - is much- 
much smaller than the thickness of the shocked plasma, most of the 
synchrotron loss occurs in the region of low magnetic field (x 2> 
Lp), which also controls the maximum Lorenz factor of electrons, 
provided that 

Bl [1 - (Lp/xof'^-'] 7(277 - 1) < Bo' (Rl/Lp), (6) 

where xo/Lp = {B^, / BqY^^ \ the above equation is obtained by 
calculating energy loss in the region of low magnetic field (Bo) 
where the electron travels a distance ~ Rl, and the loss in the 
region Lp<x < xo, and requiring the former to be larger Since 
a^o /Lp 3> 1 for the case of interest, the above condition simplifies 
to 



{Bu,/Bof <Rl/Lp 



(7) 



when rj > 1/2. Thus, one of the requirements for exceeding the 
~ 50 MeV limit is that the width of the region occupied by high 
magnetic fields (xo) is much smaller than Rl- The above condi- 
tion (eq. list also guarantees that the Larmor radius of electron in 
the region of high magnetic field is much larger than Lp, and thus 
the deflection of electron orbit while passing through this region is 
small; electrons are turned around in the region of low magnetic 
field which occupies most of the shocked plasma volume. 

The maximum LF of an electron, 'ymax, is obtained by re- 
quiring that the energy loss due to synchrotron radiation in a Lar- 
mor time (Rl/c) not exceed me^rnaxC^ /2. The case of interest is 
where electrons lose their energy while traveling through the re- 
gion of weak magnetic field. In this case ^max is same as given in 
equation l[T), i.e. 
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This Lorentz factor must satisfy the condition that Rhi'^max) is 
smaller than the comoving radial width of the shocked fluid (the 
width is R/T by causality argument), i.e. the electron is confined 
to the system. This requires 
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where F is the Lorentz factor of the shocked fluid, R is the dis- 
tance of the shock front from the center of explosion, and we have 
adopted the widely used convenient notation x„ = x/10" . 

Electrons also suffer IC loss of energy which is ignored here; 



considering that 7 

71 



10^ B, 
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IC scatterings are in Klein- 



Nishina regime and thus the IC loss is highly suppressed (Bamiol 
Duran & Kumar, 2011). 

The maximum synchrotron frequency, in shock comoving 
frame, is given by 
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(10) 

which is larger than the case of a uniform magnetic field by a factor 
/ Bq . 

We next estimate the synchrotron flux at Vmax to see whether 
it lies on the powerlaw extension of flux at lower frequencies (< 50 
MeV) or not. 

Let us consider electrons of LF > 7maa;/2^''^ which produce 
synchrotron photons of frequency > Vmaxj'i-- By the definition of 
"imax these electrons will lose half of their energy while traveling 
down-stream of the shock front. The energy fraction of these elec- 
trons lost to synchrotron radiation (^) while they travel through the 
region where the magnetic field strength is ~ B^ is 



B^RLi-Jr, 



2x10" 



Bu 



Bo 



B'^^^iFs/n, 
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where RLifmax) = rnc^^rnax/iqBo) is the Larmor radius, and 
we made use of equations Q & ([Sj to substitute for Lp & ^max- 
The frequency of radiation produced in the region of high magnetic 
field (Bm) is ~ Vmax which is given by equation dlOb . Much of the 
rest of the electron energy is lost to radiation in the region occupied 
by the weaker field Bo, and the synchrotron frequency of the emer- 
gent radiation from this region is VmaxiBo/B^) = uiow\ where 
fiom ~ 50 MeV is the maximum synchrotron frequency when the 
magnetic field is uniform. 

The specific flux at Vmax is obtained from the equation 
i^maxfii^Tnax) ~ ^i^iow f {i^iow) ■ Therefore, the spectral index be- 
tween uiow and Umax — defined as (x — is 
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ln(B„/Bo) 
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(12) 

The second equality is obtained by substituting for ^ from equation 
( II It . It can be shown, after some algebra, that jS cannot be larger 
than — p/2, where p is the index for electron distribution. 

We note that the synchrotron radiation formula can only 
be used provided that the magnetic field is coherent on scale 



Rhhry 



jrieC 



l.TxlO-" R"^cm. The co- 



herence scale for magnetic field close to the shock-front is ~ 
Lp>2xl0^(r,/ne)'/^ cm. Thus, we see that the coherence length 
for magnetic fields is much larger than RL/"/max as long as B^ ^ 
0.1 mG. The coherence scale grows with increasing distance from 
the shock front (since smaller scale fields have time to decay), and 
so does the Larmor radius. However, as long as Bo -5^ 1 mG, the 
magnetic field is coherent on length scales that are much larger than 
RL/'^max everywhere down-stream of the shock front. Therefore, 
it is safe to use the synchrotron radiation results as we have in this 
sectior0. 

When magnetic field changes on a length scale that is short compared 
with electron Larmor radius divided by its Lorentz factor, the radiation pro- 
duced is in the jitter regime which is discussed in Medvedev, (2000). 



3 APPLICATIONS TO GRBS AND SUPERNOVA 
REMNANTS 

We consider two particular applications of the results found in the 
previous section, one of which is the application to GRB afterglow 
radiation, and the other is the application to synchrotron radiation 
from SN remnants. 



3.1 Gamma-ray Bursts 

The relativistic jet produced in GRBs drives a shock wave into the 
medium in the vicinity of the burst, and synchrotron radiation from 
shock heated plasma is responsible for the long lived afterglow ra- 
diation from radio to 7-ray frequencies (eg. Piran, 2004; Meszaros 
2006; Gehrels et al. 2009; Zhang 2007). Roughly half of the GRB 
jet energy is imparted to the surrounding medium at the decelera- 
tion radius, Rd, which is given by, eg. Piran (2004), Zhang (2007) 
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(13) 

where E is the isotropic equivalent of energy carried by the jet, tiq 
is the average number density of protons in the unshocked medium 
within the region of radius Rd, and Fq is the initial Lorentz factor 
of the jet. The Lorentz factor of the jet decreases with distance as 
comoving radial width of the shocked plasma is 5R ~ 
R/r, and therefore SR/Lp ~ 10^ RnT^^ {ue/ry^^ . 

The minimum magnetic field expected in the shocked fluid is 
Bo = ATBism which is simply the field of the GRB circum-stellar- 
medium (CSM) compressed by the relativistic shock by a factor 4F; 
where Bism is the field in the CSM at a distance Rd from the 
center of explosion and is expected to be of order a few /uG. Thus, 
the condition that electrons of Lorentz factor 'jmax are confined to 
the shocked fluid by the weak magnetic field (see eq.[9]( is satisfied 
for GRB external shocks as long as the shock is highly relativistic 
(F>10). 

The field is likely amplified by a large factor immediately 
down-stream of the shock front by the two-stream instability, such 
that the energy density in the magnetic field is a fraction, es, of the 
thermal energy density of the shocked fluid (Medvedev & Loeb, 
1999). In this case 



Btu ~ (327resnompC F 



2p2Nl/2 
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with es ~ 10"^ (Medvedev & Loeb, 1999); we have used here 
the relativistic self-similar shock solution (Blandford & McKee, 
1976) according to which the thermal energy density behind the 
shock front is Anomp(?V^ . This magnetic field has a coherence 
length of order the plasma scale and decays down-stream of the 
shock front (Gruzinov, 2001; Sironi & Spitkovsky, 2011; but see 
also Medvedev et al. 2005). It is unclear what is the exact value of 
the field very far from the shock front. According to recent numeri- 
cal simulations of Sironi & Spitkovsky (201 1) the far field might be 
of order the shocked compressed ISM field; empirically we know 
from GRB afterglows that es < 10~^ for a large fraction of bursts 
(see §2. 1 for references to a number of papers on this topic). There- 
fore, for GRB afterglows we expect the ratio between B^ and Bo 
to be larger than ~ 10^ and possibly of order 



B„ (27resnomp) ' c 



Bo 



Bis 



9.5x10' (esno)'/'B-i,_5, (15) 



where Bism.,-5 is the ISM magnetic field in units of 10 ^G; the 
precise value of B^/Bo is not important as long as it is larger than 
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Figure 1. Numerical spectra for different values of Rj^/Lp and rj. The 
solid, dashed and dotted lines correspond to Rij/Lp = 10^, 3 X 10® 
and lO'', respectively (and rj = 0.7). Similarly, Ri^/Lp = 10^ and 
rj = 0.6, 0.7 and 0.9 for the red, black and blue lines, respectively. For 
all these cases, 7max/7min = 10^ ™d B^,/Bo = 10'^, and the size of 
the system is larger than ijj^ by a factor of ~ 10. The spectrum above 
Vlow = 9meC^/(167r(7^) ^ 50 MeV, has a power-law slope roughly 
equal to the spectrum at lower frequencies and it extends to well beyond 
Ulow The flux at ui^^ can exhibit a sharp decline by almost an order of 
magnitude for certain parameters, however, the flux at > IGeV even in the 
worst case lies on a slightly steeper power-law extension of the lower energy 
( < 50 MeV) spectrum. 



GeV for synchrotron radiation from shock heated plasma, eg. Piran 
& Nakar (2010). However, the >10^ MeV data for these bursts is 
in excellent agreement with synchrotron radiation from the exter- 
nal shock, which also provides a very nice fit to the late time x-ray, 
optical and radio data using the same exact parameters that provide 
a fit to the early >10^ MeV data (Kumar & Barniol Duran, 2009, 
2010). 

This conflict between the excellent fit to the high energy data 
on one hand and the conventional upper limit on synchrotron pho- 
ton energy of ~ 5 GeV on the other hand can be resolved provided 
that the magnetic field just behind the shock front is much larger 
than far down-stream; photons of energy > 5 GeV are produced 
by electrons of LF ~ ^/max while they travel through regions of 
stronger magnetic field close to the shock front. 

We note that there is a large uncertainty regarding the length 
scale over which fields decay, & Bo. However, the comoving 
frame 50 MeV limit on synchrotron photons is violated whenever 
Bn,/Bo > 10, and the distance electrons of LF "/max travel before 
losing their energy to radiation is large compared to the length- 
scale for the decay of magnetic field. The data for GRBs with > 10 
GeV emission show that the volume averaged magnetic field in the 
shocked fluid is a factor 10^ smaller than the equipartition valu^fl 
(Kumar & Barniol Duran, 2010) whereas the Weibel generated field 
strength in the immediate vicinity of the shock front is expected to 
be ~ 30% of the equipartition value (Medvedev & Loeb, 1999), 
and therefore B^j/Bo > 10^ for these bursts. 



~ 10, and the field decays on a length-scale that is short com- 
pared with the Larmor radius (Rl) or the distance the highest en- 
ergy electrons travel before losing half their energy. Substituting 
this into equation dlOl l we find 

^0.5TeV (esno)^''^Brs™.-5- 
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(16) 



The maximum energy for synchrotron photons in the observer 
frame is Tvmax/ ( 1 + . If the time it takes for an electron to travel 
from one side of the shock front to another is larger than i?L/c by 
a factor 5-10 then Umax would be smaller than given by the above 
equation by the same factor. 

Using equations il2l & l ll5t one can obtain an approximate 
analytical power-law spectral index for i^iow < i' < Vmax', how- 
ever, we have calculated the spectrum numerically for different 
choices of parameters and present the results in Figure [T] 

The spectrum in the observer frame should extend to ~ TeV 
energies provided that IC losses don't prevent electrons from accel- 
erating to ^max ~ 10^ (eq. [8]l, and these high energy photons are 
not converted to pairs while traveling from the emission site to ufl; 
the effect of IC loss on '^max, in the context of GRBs, is discussed 
in Barniol Duran & Kumar (2011) and Piran & Nakar (2010). For 
certain range of parameters the flux at photon energies much larger 
than 50 MeVxF lies below the powerlaw extrapolation of flux be- 
low viow (Fig. 1), and in those cases the flux at very high energies 
is likely to be lower than the sensitivity limit of current detectors. 

Several bursts detected by the Fermi satellite show emission 
above 10 GeV in the GRB host galaxy rest frame: Abdo et al. 
(2009a,b, 2010), Ackermann et al. (2010, 2011). These bursts vi- 
olate the conventional upper limit of ~ 50r/(l + z) MeV ~ 5 

* Photons of energy > TeV from a source at cosmological distance are con- 
verted to electron-position pairs due to interaction with infrared background 
radiation. 



3.2 Supernova remnants 

The deceleration radius for supernova remnants is 
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5xl0''*(£;5i/no)'''^U0,9''^cm, (17) 



where E is the kinetic energy of the remnant, and iio,g is its ini- 
tial speed in unit of 10^ cm/s. The radius increases with time as 



2/5 



during the adiabatic expansion phase of the remnant. The 
radial width of the remnant {5R) is roughly 1/5 of its radius, 



-1/2 



and the proton-plasma length Lp ~ 2x10 n, 
SR/Lp> 10^\ 

The equipartition magnetic field for the remnant is 



cm. Therefore, 



Bu 



(407rmpno)^''^ii 



ImG nl/'^vs, 



(18) 



where v is the remnant speed. The mean magnetic field of the 
remnant should be larger than, or equal to, the shock compressed 
magnetic field of the CSM, i.e. Bo SB^sm ~ 10/iG. Therefore, 

B»/Bo<10^ 

The maximum Lorentz factor of electron for this mean mag- 
netic field is -jmax ~ [37rg/(aTBo)]^''^ ~ 3xl0"' (see eq.[8), and 
its Larmor radius is 5xl0^*cm which is smaller than the shell width 
when the velocity drops below 10® cm/s. 

The synchrotron cooling time for electrons with Lorentz factor 
^max ~ 10^" traveling in a field of strength Bo, which occupies 
much of the volume of the remnant, is tc ^ Qi^mec/ {arB^^rnax)- 
In other words, the distance electrons travel before cooling down is 



^ The small value for volume-averaged magnetic field for these GRBs fol- 
lows from the fact that the synchrotron cooling frequency at the end of the 
burst must have been > 10 MeV; otherwise the ratio of flux at 150 keV 
to that at 100 MeV would be much larger than the observed value (large 
cooling frequency, > 10 MeV, implies small magnetic field). 
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Lc ^ etc ~ 8x10^ cm which is of order of the remnant width 
when the remnant velocity drops below lO^cms"^. 

Thus, Lc/Lp ~ 10", and so Lc/Lp > {B^/Bof ~ 10*. 
Therefore, the specific flux at 50 MeV (B^/Bq) ~ 5 GeV is a 
factor ^ ~ 10^ smaller than the flux at 50 MeV, which is perhaps 
too small to be of practical consequence. However, if the magnetic 
field in the remnant decays not on the plasma length scale but on 
~ 10* times the plasma scale then the flux at Vmax ~ 5 GeV would 
lie on the extension of the powerlaw spectrum at lower energies. 



4 CONCLUSIONS 

We have shown that photons of > 10 GeV energy observed from 
several GRBs by the Fermi satellite can be produced via the syn- 
chrotron process in the shock heated circum-stellar medium. The 
conventional wisdom that the maximum energy of photons in this 
situation should not exceed ~ 50r/(l + z) MeV ^ 5 GeV is 
violated because this limit is obtained by assuming a uniform mag- 
netic field down-stream of the shock front. However, when the field 
is much stronger close to the shock front and decays downstream 
(such as when magnetic fields are generated by Weibel-type insta- 
bilitielfl) the maximum photon energy is larger than this limit by a 
factor of the ratio of field just behind the shock front and far down- 
stream; photons of energy > 10^ GeV (in observer frame) can be 
generated this way in GRB early afterglows. 

The upper limit of ~ 50 MeV is difficult to violate for syn- 
chrotron radiation from supernova remnants unless the field were 
to decay on a length scale much larger than proton plasma length 
but smaller than the distance highest energy electrons travel before 
losing half their energy. 
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° Large variations in magnetic field strength across the shocked plasma 
could also arise as a result of field in the circum-stellar medium varying 
with distance from the progenitor star which is compressed and amplified 
near the shock front. Long-duration GRBs and SNe are produced when a 
massive star dies, and the medium within a few pai'secs of the progenitor 
star is carved out by its wind within the last ~ 10^ years of the collapse. 
Although little is known about the wind and its associated magnetic field in 
the last ~ lO^yeai's of a star's life, it is possible that the magnetic field in 
the wind could undergo lai'ge variation with time due to the magnetic cycle 
of the star driven by its rapid rotation and sub-surface convection. 
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